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(54) Magnetic resonance imaging apparatus 

(57) In a magnetic resonance imaging apparatus in 
which excitation pulses are applied to a restricted region 
5 of a magnetic resonance imaging magnet bore, in 
which the field is uniform, and the data samples collect- 
ed are Fourier transformed to form a volumetric image 



of the restricted region, a motor 6 continuously moves 
a patient couch 3 so that a region of interest passes 
through the region of good field, and the data samples 
collected are corrected to compensate for the motion so 
that a volumetric image is formed of greater length than 
that of the restricted region. 
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Description 

[0001] This invention relates to magnetic resonance 
(MR) imaging apparatus. 

[0002] In magnetic resonance imaging, the subject to 
be imaged is positioned in a strong magnetic field, for 
example, produced in the bore of a superconducting 
electromagnet, and the protons of hydrogen atoms (and 
of other MR active nuclei) in water and fat align parallel 
and anti-parallel to the main magnetic field, precessing 
around the direction of the field at the Larmor frequency. 
[0003] A transmit coil applies pulses of r.f. energy at 
the Larmor frequency in a direction orthogonal to the 
main field to excite precessing nuclei to resonance, 
which results in the net magnetisation of all MR active 
nuclei being flipped from the direction of the main mag- 
netic field into a direction having a transverse compo- 
nent in which it can be detected by the use of a receive 
coil. 

[0004] The received signal can be spatially encoded 
to produce two-dimensional (slice) or three-dimensional 
(slab) information about the distribution of MR active nu- 
clei and hence of water and tjssue. 
[0005] The received signal can be confined to a slice 
of the patient in the following way. 
[0006] Referring to Figures 1 and 2, a superconduct- 
ing electromagnet 1 (seen in side view in Figure 1 and 
in end view in Figure 2) has a bore 2 for receiving a pa- 
tient supported on a couch 3 which can be slid into the 
bore of the magnet from a position outside. By making 
the strength of the main magnetic field vary along the 
length of the bore (the z-axis) using magnetic field gra- 
dient coils to increase and/or decrease the main mag- 
netic field, it is possible to excite MR active nuclei (the 
remainder of this description refers hydrogen nuclei as 
an example) confined to a slice. orthogonal to the z-di- 
rection, since the Larmor frequency depends upon the 
strength of the magnetic field, and the frequency of the 
r.f. pulse can be chosen to correspond to that frequency 
(Gradients could equally be set up to excite slices or- 
thogonal to the x or y axes). 

[0007] Spatial encoding of the slice can be produced 
by x and y magnetic field gradient coils which alter the 
strength, but not the direction, of the main magnetic 
field, in the x and y-directions. The frequency and phase 
information in the received signal can be analysed to 
map the distribution of the.hydrogen nuclei in the plane 
of the slice. 

[0008] Slab (or three-dimensional) imaging can be 
performed by commonly exciting a region with the r.f. 
excitation pulses that is to be divided into a contiguous 
series of slices using the phase of the received signal 
to distinguish between these different slices. 
[0009] It is also known to image contiguous slabs, par- 
ticularly when attempting to trace vessels such as blood 
vessels. In one known MRI apparatus, in order to be 
able to image along the length of the leg, three separate 
slab images are recorded. The patient is moved to three 



separate axial positions for the data to be captured. The 
three images are then joined together. However, there 
are discontinuities in the images when viewed in a pro- 
jection normal to the slab boundaries, which is just the 

5 direction in which vessels would be viewed in. 

[001 0] It might be wondered why the three multi-slab 
images could not be acquired without moving the pa- 
tient, since the length of the bore of the superconducting 
electromagnet is typically much greater than the length 

io of a leg. However, the region of good field, over which 
the field is sufficiently uniform to obtain acceptably un- 
, distorted images is much smaller than the overall length 
of the bore. Typically, the. region of good field would be 
ah approximately spherical volume 4 in the centre of the 

is bore. 

[0011] While it is medically desirable to be able to col- 
lect three-dimensional image data from an elongated re- 
gion, the trend in terms of electromagnets makes this 
more difficult to achieve, since patients find shorter mag- 

20 nets to be more acceptable, and shorter magnets have 
a region of good field 5 which is compressed in the axial 
direction, although not in directions at right angles. 
[0012] The invention provides magnetic resonance 
imaging apparatus comprising means for applying r.f. 

25 excitation pulses to, and collecting volumetric data sam- 
ples from, a restricted region of the main field of the ap- 
paratus, transforming the data samples to form a volu- 
metric image of the restricted region, including means 
for advancing a patient support continuously through the 

30 restricted region, and means for correcting the data 
samples to compensate for the motion so that the volu- 
metric image formed is of greater length than that of the 
restricted region. 

[001 3] The invention enables a shorter electromagnet 
35 with a restricted length of good field to produce extended 
three-dimensional images. 

[0014] Advantageously, volumetric data samples are 
encoded with secondary phase encoding in the direction 
of motion of the patient support and with primary phase 
encoding which may have a greater number of incre- 
ments in a transverse direction, the correcting means 
being arranged to compensate successive sets of sam- 
ples of all secondary phases produced at each primary 
phase for the motion. 

[0015] Ways of carrying out the invention will now be 
described in greater detail, by way of example, with ref- 
erence to the accompanying drawings, in which: 

Figure 1 is a simplified side view of a known mag- 
netic resonance imaging apparatus; 

Figure 2 is a simplified end view of a known mag- 
netic resonance imaging apparatus; 

Figure 3 is a simplified side view of a magnetic res- 
onance imaging apparatus according to the inven- 
tion; 
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Figure 4 is a perspective view of a slab to be imaged 
in the magnetic resonance imaging apparatus of 
Figure 3; 

Figure 5a illustrates schematically the r.f. pulse of 
a pulse sequence of the magnetic resonance imag- 
ing apparatus of Figure 3; 

Figure 5b illustrates schematically the z-gradient 
pulse of the pulse . sequence of the magnetic reso- 
nance imaging apparatus of Figure 3; 

Figure 5c illustrates schematically the x-gradient 
pulse of the pulse sequence of the magnetic reso- 
nance imaging apparatus of Figure 3; 

Figure 5d illustrates schematically the y-gradient 
pulse of the pulse sequence of the magnetic reso- 
nance imaging apparatus of Figure 3; 

Figure 6a is a perspective view of a three-dimen- 
sional array of data samples in k-space; 

Figure 6b is a front view of a two-dimensional array 
of data sample in k-space corresponding to a single 
slice of the three-dimensional array of Figure 6a; 

Figure 7a shows a square object seen in plan view 
about to pass through the region of good field of the 
magnetic resonance imaging apparatus of Figure 3; 

Figure 7b shows a representation of the object in 
the x,z plane in k-space; 

Figure 7c shows a representation of the object in 
the x,z plane in the hybrid dimension of k-space and 
real space; 

Figure 7d shows the representation of Figure 7c but 
with distortion to the representation in the z-direc- 
tion for the purposes of explanation; 

Figure 8 is a plan view of slabs of a data collection 
volume passing through the region of good field in 
Figure 3; and 

Figure 9 is a block diagram showing the architecture 
of the data processing circuitry of the imaging ap- 
paratus of Figure 3. 

[0016] Like parts have been given like reference nu- 
merals throughout all the figures. 
[001 7] Referring to Figure 3, the magnetic resonance 
imaging apparatus of the invention comprises an elec- 
tromagnet 1 which is superconducting, but could if de- 
sired be resistive. The magnet has a bore 2. Those parts 
of a patient on a couch 3 which are desired to be imaged 
are translated through a region of good field of the mag- 



net 5, in which the magnetic field is sufficiently uniform 
for imaging to be possible with acceptable spatial dis- 
tortion. 

[0018] In accordance with the invention, a motor 6 is 

s provided for moving the couch 3 from a reception posi- 
tion outside the bore into the bore such that the parts of 
the patient to be imaged are translated through the re- 
gion of good field 5. And, referring to Figure 9, correction 
means 7 are provided for correcting the data samples 

10 collected to compensate for this motion. 

[0019] The superconducting magnet of the imaging 
apparatus has a main field in the z-direction and is pro- 
vided with a z-gradient coil for setting up a magnetic field 
gradient along the z-direction. X and y-gradient coils are 

is also provided to enable magnetic field gradients to be 
set up in the x and y-directions as required. 
[0020] An r.f. transmit coil 9 is provided to apply r.f. 
excitation pulses to restricted regions within the region 
of good field 5, and an r.f. receive coil 10 is provided to 

20 detect the resulting relaxation signals, which are digi- 
tised in A-D converter 11 and fed to digital signal proc- 
essor 1 2 which assembles a notional three-dimensional 
array of data in k-space. Three-dimensional Fourier 
transformation takes place in digital signal processor 1 3 

25 to convert the k-space array to an array in real space, 
and this data is passed to display 1 4. 
[0021] Referring to Figure 4, consider a slab of vol- 
ume within the region of good field 5 which it is desired 
to image. The pulse sequence is shown in Figures 5a 

30 to 5d. In each repetition time (TR), an r.f. pulse, typically 
of flip angle between 1 0 and 50° is applied to the patient 
of which the region shown in Figure 4 defines a slab, 
and the resulting echo signals are received while a mag^ 
netic field gradient in the y-frequency encode direction 

35 js being applied. This means that the echo signals, 
which consist of a combination of frequencies, can be 
analysed by Fourier transformation to detect the distri- 
bution of hydrogen nuclei in the y-direction. Accordingly, 
the y-direction is frequency encoded (Fy) in Figure 4. 

^0 The readout gradient Gy is preceded by a negative lobe 
to cancel the dephasing of spins produced during the 
readout portion. The other transverse direction, the x- 
direction, is decoded by means of phase encoding, de- 
noted by Px. Typically there may be 256 increments of 

45 phase, and a separate readout Gy is performed for each 
of them. 

[0022] The longitudinal position of the particular slice 
being imaged is determined by the z-gradient (Figure 
5b) but, since this is slab imaging, in fact a contiguous 

50 series of slices are encoded, and so the Gz gradient en- 
ergises hydrogen nuclei in the slab illustrated in Figure 
4. To distinguish one slice in the slab from the next, a 
series of secondary phase encode gradients are ap- 
plied, for each primary phase encode increment in the 

55 x-direction. That is, the slab excitation pulse Gz is fol- 
lowed by a second gradient pulse in the z-direction 
which is turned on and then off before read-out takes 
place, the gradient stepping from a negative value 
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through several Increments to a positive value. Typical- 
ly, there would be 32 such increments. Typically, Tr 
would be around 10ms. Thus, for one x increment, the 
time required would be around 320ms to run through the 
slice select phase encoding gradients. The entire vol- 
ume would be run through in a period 256 times as long, 
to allow each of the x phase increments to be executed. 
In other words, the whole volume shown in Figure 4 
would be scanned on this basis in around 80 seconds. 
[0023] The time varying echo samples do not directly 
represent hydrogen nuclei distribution, but fill a three- 
dimensional array in so-called k-space shown in Figure 
6a, where dimension Ky represents frequency and di- 
mension Kx and Kz represent frequency synthesised 
from phase encoding. Referring to Figure 6b, which 
shows one of the slices of the array 6a, each echo signal 
produces a series of points on one vertical line Ky, and 
a fresh vertical line is produced at each Kz phase en- 
code increment. 

[0024] The k^space of Figure 6a is converted to a rep- 
resentation of the distribution of hydrogen nuclei in real 
space in digital signal processor 1 3, by separable Fou- 
rier transformation along each of the x,y and z axis. 
[0025] The relationship between real space and k- 
space can be further explained by reference to Figure 
7a to 7d. Thus, if one were to imagine that the object 
being imaged was a square-sided figure, the sides of 
which are aligned with the x and z-axis in Figure 3, in 
fact the k-space representation of this would appear as 
in Figure 7b. After Fourier transformation along the z- 
axis, the representation in hybrid Kx, real z-space would 
appear as in Figure 7c, where the array would now con- 
sist of points arranged on lines parallel to the z-axis and 
equal in number to the number of k-phase encode in- 
crements. Fourier transformation along the x-axis would 
yield the square of Figure 7a. 

[0026] This correspondence between real space and 
k-space, and Fourier transformation back to real space 
is of course well known. However, what is not known is 
for the patient to be translated through the region of 
good field 5, in the direction of the axis z, during the proc- 
ess of data collection. 

[002*7] The effect of this can be considered by refer- 
ence to Figure 7a. Let us image that the time taken to 
step through all the phase encode increments in the z- 
direction is sufficiently small that the motion of the four- 
sided figure can be ignored. Then, if z^ and z 2 are the z 
co-ordinates of the top and bottom limbs of the figure at 
the first x phase encode increment, say, the most neg- 
ative, then the top and bottom limbs of the figure will 
have moved to the position z 3 , z 4 by the time data is 
collected at the most positive x phase encode incre- 
ment. 

[0028] While the data will still form a three-dimension- 
al matrix with orthogonal axes and will therefore follow 
the pattern of Figure 7c, one could distort the represen- 
tation of Figure 7c into that of Figure 7d in order to illus- 
trate the translation which has taken place in the four- 



sided figure between the most negative primary phase 
encode increment (z 2 ) and the most positive (z 3 ). 
[0029] If one considers Figure 7c for the case when 
the four-sided figure is stationary, Fourier transforma- 

s tion in the x-direction to reconstitute the right hand limb 
of the four-sided figure requires computation using each 
point at the right hand end of each of the Kx lines. In 
view of the movement which has taken place while data 
is being collected, the points which now need to be com- 

10 bined are those under the slant line z 2 ,z 3 in Figure 7c. 
Knowledge of the speed of translation of the couch en- 
ables this spatial correction to be readily calculated, and 
this is done in correction means 7. The result of this is 
that the data in the hybrid space of Figure 7c can be 

is corrected for the motion of the couch. 

[0030] In fact, the point z 3 in Figure 7c, as well as rep- 
resenting a displacement related to the speed of move- 
ment of the couch, also reflects a phase change de- 
pendent upon the speed of the couch. It is then merely 

20 necessary, using the Fourier shift theorem, to multiply 
each point in the k-space array of Figure 6a by an ap- 
propriate phase correction to take account of the trans- 
lation of the volume while the data was being collected. 
If this phase correction is made, then the correction can 

25 be performed wholly in k-space. The spatial correction 
discussed with relation to Figure 7c on the other hand 
requires that the z-dimension is first of all Fourier trans- 
formed to real space. A description of the use of the Fou- 
rier shift theorem used for correction of involuntary mo- 

30 tion is magnetic resonance imaging is described in Fre- 
quency Domain Simulation of MR Tagging, William R 
Crum et al, J MR 1.1 998; 8:1040-1050. 
[0031] Considering the three-dimensional slabpf Fig- 
ure 4, if the x,z-dimensions coincide with those of the 

35 four-sided figure shown in Figure 7a, if the gradient Gz 
were arranged to excite the entire volume between z-, 
and z 2 , it will be apparent that only the region of the slab 
between the coordinates z 1 and z 3 would be fully imaged 
(all primary and secondary phase encodes acquired), 

40 since the region of the slab lying between z-coordinates 
z 4 and z^ would be outside the excitation volume (be- 
tween z n and z 2 ) at the end of primary phase encoding. 
[0032] In reality, the volume to be imaged will extend 
further in the z-direction than the slab of Figure 4. Thus, 

45 after the first scan has been completed, the volume of 
the patient originally lying between z-coordinates z 1 and 
z 2 has been translated to z-coordinates z 4 , z 3 . The vol- 
ume of the patient lying between z-coordinates z-, and 
z 3 has been imaged, and a fresh volume of the patient 

50 between z-coordinates z 3 and z 2 has now entered the 
volume which is excited (between z^ and z 2 ). Conse- 
quently, after the next scan has been completed, the 
fresh volume between z 3 and z 2 is imaged, and so on. 
[0033] Of course, with the coordinates shown in Fig- 

55 ure 7a, a large proportion of the patient's volume which 
had been imaged on one scan would be imaged also on 
the next scan. Duplicated data would be averaged to 
improve the s:n ratio. Accordingly, the excitation volume 
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and the speed of translation may be chosen so that point 
z 3 moves all the way to point in the time which one 
scan takes to perform. The volume which is fully imaged 
on each scan is now just a single slice, and there is no, 
overlap between the volume fully imaged on one scan s 
and the volume fully imaged on the next scan. Continu- 
ous imaging of a patient on the couch 3 moving through 
the region of good field takes place, and the result is a 
continuous seamless length of time 3D image data. 
[0034] It turns out that a certain amount of overlap be- 10 
tween successively imaged volumes can in fact be de- 
sirable. Thus, the scheme is most sensitive to errors 
where the join in the realigned data in hybrid Zk* space 
is close to the central K x =0. To guard against artifacts 
produced by errors in this region it is helpful to use data 15 
interpolation to ensure a smooth join or to double collect 
/>data so that the central data always comes from a 
narrow range of positions within the slab. This can be 
achieved by making the translation velocity smaller so 
that a distance less than Z-, to Z 3 is traversed in the time 20 
taken to acquire all primary and secondary phase en- 
code steps, or by acquiring the central lines twice, for 
example if sequentially ordered phase encoding is used 
the central lines would he acquired in the natural order 
and again at the end of the series. As another alterna- 25 
tive, referring to Figure 8, which is a plan view of four 
slabs A, B, C, D, successive frequency offsets could be 
applied to the frequency of the r.f. pulse, so that in the 
first data collection, the data in volume A was collected. 
The frequency of the r.f. pulse could then be offset to 30 
collect data in volume B. It could then be offset of collect 
data in volume C. It could then be offset to collect data 
in volume D. The r.f. frequency could then be reset to 
collect data in volume A again etc. 

[0035] It is possible that there would be discontinuities 35 
between the slabs produced in this method. This could 
be reduced by "walking" a thin slab along the z-axis con- 
tinuously while acquiring data in an interleaved fashion 
along the Ky axis (Sliding Interleaved ky (SLINKY) Ac- 
quisition: A Novel 3D MRA Technique with Suppressed 40 
Boundary Artifact, JMRI 1998, Volume 8, No. 4, pages 
903 to 911, and European Patent Application No. 
98308921.0). 

[0036] The above discussion has taken place with the 
assumption that the motion of the object being imaged 45 
is insignificant over one run through of all the secondary 
phase encode gradients. The assumption holds provid- 
ed that the sloping line in Figure 7c does not cross the 
second from the bottom primary phase encode line be- 
yond the second data point of that line. However, the so 
effects of the couch motion while a single run through 
of secondary phase encode data takes place can in fact 
be corrected for. in the same way that the motion be- 
tween the lines was visualised in Figure 7c, the motion 
during a line could be visualised in Figure 7c by slightly 55 
displacing each of the points forming a z line. In other 
words, a velocity correction could be applied to each 
point, for example, in the form of a phase correction. 



[0037] It has been assumed that a separate r.f. pulse 
is required for each secondary phase encode incre- 
ment. Of course, it would be possible for a multi-echo 
sequence to be used, so that a different secondary 
phase encode gradient was applied before each repeat- 
ed echo. Suitable pulse trains would be of the generic 
form termed multi-shot echo planar imaging with a sec- 
ondary phase encode or multi-shot echo volumar imag- 
ing. 

[0038] The description has referred to data being as- 
sembled in orthogonal arrays in what would be cartesian 
coordinates, but the manipulation of the data could 
equally take place in polar coordinates. 
[0039] The invention is particularly applicable to mag- 
nets having only a narrow region of good field in. the z- 
direction eg. maximum thickness of 3cms. This would 
apply particularly to shorter magnets. While an annular 
magnet has been described, the invention is also appli- 
cable to other geometries for example C-magnets and 
including permanent magnets. 

[0040] The invention is of course applicable to slice 
select and phase encode gradients being in the x or y 
or oblique directions, with appropriate changes to the 
motion of the volume to be imaged. This invention can 
also include correction for spatial distortions introduced 
by main field inhomogeneity and/or gradient non-linear- 
ity. The correction could take the form of a real space 
unwarping or a k space-regridding. 



Claims 

1 . Magnetic resonance imaging apparatus comprising 
means for applying r.f. excitation pulses to, and col- 
lecting volumetric data samples from, a restricted 
region of the main field of the apparatus, transform- 
ing the data samples to form a volumetric image of 
the restricted region, including means for advancing 
a patient support continuously through the restrict- 
ed region, and means for correcting the data sam- 
ples to compensate for the motion so that the volu- 
metric image formed is of greater length than that 
of the restricted region. 

2. Magnetic resonance imaging apparatus as claimed 
in Claim 1, in which volumetric data samples are 
encoded with secondary phase encoding in the di- 
rection of motion of the patient support and with pri- 
mary phase encoding in a transverse direction, the 
correcting means being arranged to compensate 
successive sets of samples of all secondary phases 
produced at each primary phase for the motion. 

3. Magnetic resonance imaging apparatus as claimed 
in Claim 2, in which the correction means is ar- 
ranged to apply a motion correction to rows of data 
samples in hybrid real space/k-space. 



5 



9 EP1 024 371 A2 

4. Magnetic resonance imaging apparatus as claimed 
in Claim 2, in which the correction means is ar- 
ranged to apply a phase correction to rows of data 
samples in k-space. 

5. Magnetic resonance imaging apparatus as claimed 
in any one of Claims 2 to 4, in which the primary 
phase encoding has a greater number of incre- 
ments than the secondary phase encoding. 

6. Magnetic resonance imaging apparatus as claimed 
in Claim 1, in which volumetric data samples are 
encoded with secondary phase encoding in the di- 
rection of motion of the patient support and with pri- 
mary phase encoding in a transverse direction, the 
correcting means being arranged to compensate all 
samples individually. 



20 



25 



30 



35 



40 



45 



50 



10 



6 



EP 1 024 371 A2 




7 



EP 1 024 371 A2 



R.F. 



i 



FIG. 5a 



SLICE SELECT 
AND SECONDARY 
PHASE ENCODE 




FIG. 5b 



n - 






* PRIMARY PHASE 








ENCODE 









FIG. 5c 



y FREQUENCY ENCODE 
(READ-OUT) 




FIG. 5d 




FIG. 6a k 



FIG. 6b 




FIG. 8 



8 



EP 1 024 371 A2 



REAL SPACE 



K - SPACE 




FIG. 7a 




FIG. 7b 




FIG. 7c 



FIG. 7d 



9 

2. 



8 



RF 

TRANSMIT 
COIL 



GRADIENT 
COILS 



10 



T 
11 



RF 












RECEIVE 




A/ D 




DSP 




COIL 













6 



MOTOR 

ZD 



CORRECTION 
MEANS 



7 V 

12 FIG. 9 7 



- DSP 



S 

13 



- DISPLAY 



14 



9 



